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Abstract: 

A new method for the reduction of one-loop tensor 5-point integrals to related 4- 
point integrals is proposed. In contrast to the usual Passarino-Veltman reduction and 
other methods used in the literature, this reduction avoids the occurrence of inverse 
Gram determinants, which potentially cause severe numerical instabilities in practical 
calculations. Explicit results for the 5-point tensor coefficients are presented up to rank 4. 
The expressions for the reduction of the relevant 1-, 2-, 3-, and 4-point tensor coefficients 
to scalar integrals are also included; apart from these standard integrals no other integrals 
are needed. 
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1 Introduction 



High-energy collider experiments reached a high level of accuracy in the last decades. 
For instance, the LEP, SLC, and Tevatron experiments tested the Standard Model of 
electroweak and strong interactions as quantum field theories, in the sense that quantum 
corrections, i.e., higher-order perturbative radiative corrections, had to be taken into 
account to successfully compare predictions with data. The most precise results typically 
were obtained from investigations of 1 — > 2 particle decays and 2^2 scattering reactions. 
At present and future colliders, such as the Tevatron, the LHC, and an e + e~ linear collider, 
scattering reactions with more final-state particles will gain increasing interest in various 
contexts, such as in the production of jets, heavy quark flavours, electroweak gauge bosons, 
and Higgs bosons. 

In this paper we focus on 2 — > 3 particle reactions and present one of the basic 
ingredients in the evaluation of radiative corrections at the one-loop level, the calculation 
of 5-point tensor integrals. Following the well-known procedure of Passarino and Veltman 
[1], one-loop tensor integrals can be recursively reduced to scalar integrals by solving sets of 
linear equations for the tensor coefficients in the most general Lorentz-covariant ansatz. 
Specifically, for each tensor rank this procedure involves a factor of the inverse Gram 
matrix which is built by the momenta that span the tensor. Since these momenta become 
linearly dependent at the boundary of phase space, the inverse Gram matrix becomes 
singular at the boundary. In practice, this leads to numerical instabilities near the phase- 
space boundary, where the zero in the Gram determinant appearing in the denominator 
of a tensor coefficient is compensated by a delicate numerical cancellation between scalar 
integrals in the numerator. The situation is not very problematic in 2 — > 2 reactions, 
where the boundary is reached in forward and backward scattering, i.e., in isolated points 
in a single phase-space variable (scattering angle 6 — 0, ir). In the three-particle phase 
space the situation is more involved, and the numerical problems caused by vanishing 
Gram determinants of four momenta turn out to be serious. Fortunately the appearance 
of these Gram determinants can be avoided completely if the four-dimensionality of space 
time is exploited. In the following we describe such a procedure. 

Other methods for calculating tensor 5-point integrals via reduction to simpler inte- 
grals have been described in Ref. [2]; however, in these approaches either Gram determi- 
nants were not avoided completely, or the set of standard integrals had to be extended. 
Various treatments of scalar 5-point integrals have been presented in Refs. [2-6]. Methods 
for the numerical integration of scalar 5- and 6-point integrals [7] and of general one-loop 
integrals with up to six external legs [8] with special treatments of singularities have been 
proposed recently. 

It is known since a long time that the scalar 5-point integral can be reduced to scalar 
4-point integrals by using the linear dependence of the integration momentum on the four 
external momenta in four space-time dimensions [3,4]. In Ref. [5] a generalization of this 
reduction formula to tensor 5-point integrals has been proposed, but some extra terms 
were missing in the derivation. We supplement the derivation of Ref. [5] by these missing 
terms and work out explicit formulas for the tensor coefficients up to rank 4. The relevant 
formulas for the reduction of the coefficients of the 1-, 2-, 3-, and 4-point functions to 
standard scalar integrals are also presented. Apart from these standard functions no 
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integrals are needed. In summary, the results of this paper comprise all one-loop tensor 
integrals that occur in 1 — > 2, 2 — * 2, and 2 — » 3 particle reactions with up to four external 
gauge bosons (and an arbitrary number of spin-| fermions and scalars) in renormalizable 
gauge theories. Apart from the six-point functions, the presented set covers also all one- 
loop tensor integrals appearing in processes with six external particles including up to two 
external bosons, i.e., in particular, the processes 2f — > 4/ and 2f — > 2f + 2V. Only scalar 
1-, 2-, 3-, and 4-point integrals have to be calculated in addition; to this end, methods 
and explicit results can be found in the literature [2,5,9,10]. 

In its present formulation, the described approach applies to tensor 5-point integrals 
in four space-time dimensions, i.e., possible infrared (IR) singularities of soft or collinear 
origin are assumed to be regularized by off-shell or mass regulators. 1 This procedure is 
widely used in the calculation of electroweak corrections. However, the present formalism 
can also be used in dimensional regularization, which is usually adopted in QCD for the 
treatment of IR singularities. The translation of mass regulators into dimensional regular- 
ization is generally described in Ref. [11] for complete QCD and SUSY-QCD amplitudes. 
A possible strategy for the corresponding treatment of individual integrals can be found 
in Ref. [12]. 

The method proposed in this paper has already been tested in actual calculations: 
it has been used in the calculation of the QCD corrections to the process gg/qq — > ttH 
[12] and the electroweak corrections to e + e" — > i^H [13]. The results are in numerical 
agreement with an evaluation that is based on the usual Passarino-Veltman reduction, 
but a drastic improvement in the numerical stability of the 5-point function was observed. 

The paper is organized as follows. In Section 2 we derive the general relation between 
tensor 5-point integrals and the related 4-point functions. The explicit results for the 
tensor coefficients of the 5-point integrals up to rank 4 are presented in Section 3. In 
Section 4 we describe some consistency checks and applications of the presented results; 
details on the generalizations of the method to dimensional regularization can also be 
found there. Section 5 contains our conclusion. In App. A we describe the derivation of 
extra terms that arise from UV divergences in intermediate expressions. Appendices B 
and C comprise useful results for 1-, 2-, 3-, and 4-point functions. 



2 Derivation of the reduction formula for tensor 5-point integrals 

In the reduction of tensor 5-point integrals to 4-point integrals we closely follow the 
strategy and notation of Refs. [5,14]. The one-loop 4- and 5-point functions are defined 

as 

D{Q^ ilv ^ V p, ilV pa^vpaT} (Pi, Vli P3, m 0, m l, m 2, Wi 3 ) 



/dV 



ivr 2 J 1 N Q N X N 2 N Z 
E{o,p,pu,pup,^upa} (Pi,P2,P3,P4,m ,m 1 ,m2,m 3 ,m 4 ) 
{1, QnQu, q^q v q P , qpqvq P q*} 



— 2 [d A q 
m z J 



NqN^NsN^ 



(2.1) 



1 In 5-point functions, ultraviolet (UV) divergences appear only for tensor rank > 6, which is not 
considered in this paper. 
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with the denominator factors 



N i = {q + Pi) 2 - m \ + ie, 



0,...,4, Po = 0, 



(2.2) 



where ie (e > 0) denotes an infinitesimal imaginary part. Apart from D^ upa and D iiVpaT ^ 
all above tensor integrals are UV finite and can be evaluated in four dimensions (D — 4). 

The reduction of the 5-point function to 4-point functions is based on the fact that 
in four dimensions the integration momentum q depends linearly on the four external 
momenta Pi [3] . This gives rise to the identity 



2q 2 2qpi 
2p x q 2pipi 



2qp A 
2pvp A 



2p A q 2p A p x . . . 2p A p A 



2N + Y 00 2q Pl 
N 1 -N + Y w - Y 00 2p lPl 

N A -N + Y m - Y 00 2p APl 



2qp A 
2p\p A 

2p A p A 



(2.3) 



with 



Yij = m1 + m) - {pi - pj) 2 , i, j = 0, . . . , 4. 



(2.4) 



Equation (2.3) implies 



o = - 2 /A 

171 2 



9w 



-A 2 



A^oA^! • • • A" 4 q 2 - A 2 



2A^o + Y 00 2q Pl 
Ni - N + Y w - Y 00 2p lPl 



2qp A 

2pm 



N A - N + Y A0 - Fqo 2p A pi . . . 2p A p A 



(2.5) 



where P denotes the number of integration momenta in the numerator. UV divergences 
that occur in intermediate steps by expanding the determinant are regularized temporarily 
with a cutoff A — > oo, in order to be able to exploit the four-dimensionality of space-time. 
This approach is valid for all considered 5-point integrals, i.e., for P < 4. Expanding the 
determinant along the first column, we obtain 



2pipi . . . 2pip 4 



E(-ir{ D 



(A) 
a/ii. 



i=l 



X 



2p1 . 
2pipi . 

2pi-\pi . 
2p i+ ipi . 

2p A p x . 



hp' 



2p A p x . . . 2p A p A 
-D^l.., p (0) + (Y iO -Y oo )E^ p } 



2p\p A 

2pi~ip A 
2p i+ ip A 

2p A p A 



(2.6) 
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where (i) denotes the 4-point function that is obtained from the 5-point function 

Ej^ Mp by omitting the ith propagator Nf 1 . The superscript (A) indicates the regulariza- 
tion as introduced in (2.5). Since all appearing 5-point functions are UV finite (rank < 5), 
we can directly perform the limit A — > oo there and omit the superscript (A) for all E 
functions. The same can be done for the UV-finite 4-point functions. The UV-divergent 
4-point functions are for asymptotically large A split as follows, 



£>(A) _ 



£)(fin) + A 



>oo Hi— Hp 



Hi—Hp i 



is the UV-finite part of the usual 4-point function D 



where D^ u 

Hi— Hp 

defined in (2.1). The UV divergence is subtracted from D 
Whenever the functions D 



Hi -Hp 



(2.7) 



in D dimensions 



Hi— Hp 



as in the MS scheme. 2 



Hi— Hp 



(or their coefficients) are finite, D 



Hi ■ 



■hp and d{ hi-hp are 

identical, and we simply write D^ l 41p . The terms A Ml ... MP contain the dependence on 
A, but are also finite for finite A. Inserting (2.7) into (2.6) and taking A asymptotically 
large, results in 



Hi— HP 

(0) + Y 00 E i 



Hi— HP 



2pipi . . . 2p x p 4 



+E(-i) i K; 1 ) 



■HP 



D - 



i=l 



2p 4 p x . . . 2p 4 p 4 



+ A 



ani---HP 



(0 " A, 



otni—HP 



(0) + (Y l0 - Y 00 )E, 



a/ii -HP 



2p« 
2pipi 

2pj-iPi 
2p i+ ipi 

2p 4 p 1 



2pl 
2pvp 4 



2pt-m 

2p i+ lP4 

2p 4 p 4 



• (2-8) 



The terms involving Pia-D^ n .^ p (0) have been added for later convenience. The extra terms 
involving A^...^ are absent for P < 2 where no UV-singular integrals appear; they drop 
out in the final result for P = 3, but contribute for P = 4. 

The contributions to 4-point functions that involve a momentum pj a will be simplified 
in the following. To this end, we introduce the Lorentz-covariant decompositions 



n(fin) 

-^afii-.-HP 



(0 



PS... w ,(oi (p) 



4 

Pja-^j,Hi-HP (*)> 



(01 



(9) 



o, 



2 The final results for the 5-point functions E lll ,,, flp do not depend on the details of this subtraction, 
because these results are finite. We perform the subtraction in order to avoid confusion in products 
involving metric tensors. Since the subtraction guarantees that all appearing coefficients are finite, the 
question of whether we take the metric tensor in 4 or D dimensions concerns only irrelevant terms of 
0{D- 4). 
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[^ani...np(^)]^ ~ E 9ouijYj, l i 1 ...ij, j _ 1 ij lj+1 ...iAp 1 



0, 



i=2 



(2.9) 



The operation "(<?)" isolates all tensor structures in which the first Lorentz index appears 
at a metric tensor; the remaining part of the tensor furnishes the "(p)" contribution. The 
last decomposition in (2.9) becomes obvious after performing a shift q — > g— pi in the inte- 
gral. From (2.9) it follows immediately that the terms in (2.8) that involve ^ P (2)P 
drop out when multiplied with the determinants, because the resulting determinants van- 
ish. Similarly, the contribution proportional to ^(0) + Pia^S n .^ p (0)P vanishes 
after summation over i. Finally, the term Pi a Dj^}„ p (0) contributes only for i — 1, where 
it can be combined with the first term in (2.8). Rewriting the resulting equation using 
determinants and reinserting the explicit form of the tensor integrals leads to 3 



(27T/X) 



4-D 



17T Z 



/ 



d D 



Q 



Qjh ■ ■ ■ Qfj.p 
NoNi ■ ■ ■ N 4 



No 



10 



2qp 4 

2pm 



Y 40 - loo 2p 4 Pi • • • 2p 4 p 4 



(fin) 



V, 



111... tip 



+ u, 



Hi-HP- 



Here we introduced 



V, 



2p« 



2p2 



V 



afj,i...fjp 



(1) 2pipi . . . 2pip 4 



V, 



afii...fip 



(4) 2p A p 1 . . . 2p 4 p 4 



-l)^det(4 4) )2^, 



'a/ii.../i P 



(2.10) 



(2.11) 



with 



which collects all terms involving [D]^, 
the scalar 5-point integral since D a = [D c 
from the four-dimensional Gram matrix 



£)(fin) 



. .^WP- < = 1,. ..,4, (2.12) 

The extra term V is absent for the reduction of 
](p)_ The three-dimensional matrices result 



/ 2pipi . . . 2p l p 4 \ 



^ 2p 4 p x . . . 2p 4 p 4 J 



(2.13) 



by discarding the ith row and jth column. The term {7 Ml ... MP resulting from the UV 
divergences reads 



U, 



m—iip 



-2A 



..,, ; .(())<l<utZ ") + £ (-1)^' det(^ ) )2^[A am ..., p ( 



^a/ii...// P (0)]. 

(2.14) 



3 In Ref. [5] the V" and U terms are missing. 
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The actual evaluation of £/ w „ p is described in App. A and yields 



U, 



Ml— MP 



for P < 3, 

— IfK^ mim2$ M3M4 + 9\i\iiz9n±ni + 9 tinj,4 9 112 ^3) det(Z^ for P = 4. 



(2.15) 



We did not investigate these terms for P > 5. 
Using 

2piPj = Yij - Y i0 - Y 0j + Y 00 , 2q Pj = Nj - N + F 0j - loo, 



(2.16) 



we can transform the determinant on the left-hand side of (2.10) by adding the first 
column to each of the other columns, and then enlarging the determinant by one column 
and one row, resulting in 



1 Y 00 

-mp 

Y 10 — ^00 

Equation (2.17) is equivalent to 



*04 

••mp 

^14 — Y 4 



Y44 — Y ( 



04 



^Mi.-MP + ^J[X\...jlp- 



(2.17) 



^,..,p -D$i P (o) -^a,(2) -^a P (3) -^£ n U4) 



1 




Y01 


Y02 


Y)3 


Y 4 


1 


Y 10 


Y n 


Y 12 


Y 13 


Y u 


1 


Y20 


Y 2 i 


Y22 


Y23 


Y24 


1 


Y30 


Y 31 


Y32 


Y33 


Y34 


1 


Y m 


Y 41 


Y 42 


Y i3 


Y 44 



v„ 



(2.18) 



"Mi— MP ~ ^Mi— MP' 

which expresses the tensor 5-point function E^ 14lp in terms of five tensor 4-point functions 

4 A^fvW\ 



E 



>i. -MP 



_ * det(y) (fin) * i+j detjZjj , 

det(F) m - MP i,>« ; \> J ^"/'i-mp 

1 



+ 



det(y) 



(0 



det(y) 



17 



M1--MP ' 



(2.19) 



where y = (y?), i = 0, . . . 4, and y is obtained from the 5-dimensional Cayley matrix Y 
by replacing all entries in the ith column with 1. 
With the shorthand notation 



Vi 



det(y) 
det(y)' 



2 = 0,. ..,4, 



& = (-i)*"*^ 



(4), 



det(y) 



0*, i,j = l,...,4, (2.20) 
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this reads 

E^.., P = -EliO) + 2 E CijP?V a » 1 ...» P (i) + -r^U„ 1 ...„ p (2.21) 
i=o i,j=i aeiyi ) 

with U given in (2.15) for P < 4. 

3 Explicit formulas for 5-point tensor coefficients 

Here we further exploit (2.21) to derive explicit formulas for the coefficients of tensor 
5-point integrals appearing in convenient decompositions into Lorentz covariants. In order 
to be able to write down the tensor decompositions in a concise way we introduce the 
notation 

rp\)J.\...y,p\ _ rpfj, 1 ...fj,p _|_ rpp 2 ---PPPl _J_ _|_ J^PMl—MP-l ^3 ]^ 

i.e., Lorentz indices within square brackets represent a sum over all tensors with cyclic 
permutations of these indices. For example, we have 

9[pv9p\a = 9pv9pu + 9up9p<7 + QpfiQua- (3.2) 

This notation can be iterated, e.g. , 

j[/»M] _ y[/wp] _|_ y[/*H _ rppvp _|_ j>i/p// _|_ yp^^ _|_ ^wz/ _|_ yp^p _|_ j.i//ip ^ 2^ 

After cancelling the denominator N the resulting tensor integrals are not in the stan- 
dard form but can be expressed in terms of standard integrals by shifting the integration 
momentum. We choose to perform the shift q — > q — pi, so that the following 4-point 
integrals appear: 

n / n \ _ (27r/i) (4 ~ g) r D {1, g^, g^, g^g^, g^g^r} 

^{o,^ w , m ^r}W yog N^N^ 

Ni = {g + Pl - Pl ) 2 - m 2 t + ie, i = l,...,4. (3.4) 

Note that the scalar integral -D , an d the tensor coefficients -D o an d -Doooo are invariant 
under this shift. Therefore, we can omit the tilde on these functions. In the decomposition 
of Dfanv^vp^vfxj^vptrT^i) with i — 1, . . . ,4 shifted indices appear which we denote as 

ji=l j iOIl>h (3.5) 
[j -1 for % < j. 

We also use the notation 5ij = 1 — 5ij, i.e., J2i ■ •) — ■ •)• 

(i) Scalar integral 

For the scalar 5-point function (2.21) reads 

4 

£ = -5>A)(i), (3.6) 



i=0 



as was also obtained in Refs. [3-5,14]. 
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(ii) Vector integral 

For the vector 5-point function and the relevant 4-point functions we have the following 
covariant decompositions 

£>"(<) = E^(0, < = i,...,4, 

2 = 1 

D"(0) = D"(0) -p^o(0) = Efe -Pi)"^-i(0) -P^o(0), 

J'=2 

^"(i) =^[A)o«-A)o(0)]. (3.7) 
Inserting these into (2.21) for the vector integral, we find for the components 

Ej = -Y.mD^h - VoDjiO) + 2£Oi[A)o« - A)o(0)], J = 1, . . . ,4, (3.8) 
i=i i=i 

with 

A(0) = -EA'(O) - A)(0), D,(0) = A-i(O), j = 2,3,4. (3.9) 
j'=i 

(raj Tensor integral of rank 2 

For the case of the second rank tensor we introduce the covariant decompositions 

j,k=i 

D^(i) = E P»PlD jih (i) + (TA»(i), < = 1, . . • , 4, 

j,fc=l 
2i*^» 

D^(0) = L>^(0) -p^(O) +p^D o (0) 

= E (pj - - Pi)"A--i,fc-i(o) + <r Ax>(o) 

j,fc=2 

- Eri M fe -pi)^ ] A-i(o) +rfK^o(o), 

i=2 



£>°^(i) = [D a ^(i) - £^(0) 



(9) 



E^^ooi.W - E^fe -Pi)" ] ^ooj-i(0) + p 1 V Ia A»(0). (3.10) 



(9) 



i=}. i= 2 

2#» 
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The Eqq term in the decomposition of E^ v is redundant in the sense that the corresponding 
covariant, g^ u , can be expressed by the four linearly independent vectors pj up to terms 
of 0(D -4), 

ST = E ^PU zi %k + - 4). (3.11) 

j,k=i 

This means that the set of covariants in the decomposition of E^ v is overcomplete, and 
the coefficient E 00 can be defined by convenience. We use this freedom to avoid det(Z^) 
in the denominators. 

With the above decompositions we find from (2.21) 



^oo = -J2Vi D oo(i) - r] D 00 (0), 
i=i 

E jk = 2 ^ECjiiDook^hi - D 00k (0)] + (j <-> fc) J 

4 

-J2Vi D jiki(i)8jihi-VoD jk (0), j,k = l,...,A, (3.12) 



with 



i=l 



3 3 

£> u (0) = E A"*(0) + 2EA-(0) + ^o(0), 
j,fc=i j=i 

Dji(0) = -E^(0)-Oj-i(0), 
fe=l 

£> jfc (0) = 4-_i, fc _i(0), 

3 

^ooi (0) = -E Axy(O) - A)o(0), 
j'=i 

Dooj(O) = A)oj-i(O), j,k = 2,3,4. (3.13) 

(iv) Tensor integral of rank 3 

The covariant decompositions of the integrals appearing in the reduction of the 3rd 
rank tensor read 

E^= E PjPlPiE.ki + E g^PjEooj, 

j,k,l=l j=l 

m(i)= e ^rfrf^W + E^pf^W, i = i,---,4, 

j,k,l = l 2 = 1 

4 

= E (pj - - pi)"(Pi - Pi)"4--i,fc-i,t-i(o) 
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+ Ert,-Pi) p] A)o^i(o) 

" pY E (Pi - PiYiPk ~ Pi) p] Dj-i, k -M - A)o(O) 

j,fc=2 

+ pfe Efe - pi) p] A--i(o) - p»o(o), 

i=2 



j,k=i 



• E ^(P; -PiYiPk -PiY ] Axw-i,k-i(0) 

j,k=2 

4 



+ P [ i Efe - Pi) V 11q A> Oj -i(0) - s^KK 1 A)o(O). (3.14) 

i=2 

Again the coefficients -Eooj are redundant and introduced for convenience. The coefficients 
of the 3rd rank tensor 5-point function are then given by 

i=i i=i 

^yw = 1 2 Y.Cji[ D ook i i i (i)hiSu - D om (0)] + (j <-> fc) + (j <-> Z) I 
^ i=i J 

4 

-^ViDjihkiWjihiSii-VoDjkiiO) j,k,l = l,...,4, (3.15) 



with 



i=i 



3 3 3 

An(o) = - E A-«(o) - 3 E A-fc(o) - sE A(o) - ^o(o), 

j,fc,/=l j 5 fc=l j=l 

3 3 

D, n (0) = E 4"-i,h(0) + 2 E 4"-i,fc(0) + A'-i(O), 
fe,j=i fe=i 

^i(O) = -^^^(O)-^^^), 

!=1 

^■«(0) = £,-_!,*_!,,_!(()), 

3 3 

A)on(O) = E Axyk(O) + 2^ooj(0) + A)o(0), 
j,fc=i j=i 

Axyi(O) = - E A)oj-i,*(0) - Axw-i(O), 
fc=l 
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Axyfc(O) = Ax>j-i,k-i(0), j, k, I = 2,3,4. (3.16) 

(v) Tensor integral of rank 4 

For the reduction of the 4th rank tensor we introduce 

j,k,l,m=l j,k=l 

j,k,l,m—l j,k=l 

+ g^g p] °Dgll(i), i = l,..., 4, 

fl TO ( fin '(0) = £>m(fin)( ) _ p ^H(o) +p^ J D" ]a (0) +p^£>"" ] (0) - pfViK ^ (0) 
+ p1PiP P iP a iD (0) 

4 

= E (Pj - Pl)"(P* - PlTiPl - PlYiPm - Pl) <T -Dj-l,fc-l,I-l,m-l(0) 
j',fc,Z,m=2 

+ E fo^te - pi) p] fe - + 9^ - piTiPk - pi) p] ]D o,-i,k-i(o) 

j,k=2 

+ g^g p]a D^l(0) -p [ r E (Pi -PiYiPk -PiYipi -Pi) CTl 4-i,^-i(0) 
-^E^fe-pi)^ Axy-i(O) 

J'=2 

+ E [pfe(Pi -PiY\pk -piT +P a 1 P [ f(p j -PiTiPk -Pi)" ] ]A--i,fc-i(o) 

j,fc=2 

+ (piW 1ct +PiPiV p] )A) (o) 

- PiVK Efe " Pi) ffI A--i(0) + rfPifiP'MO), 

jjanvpa^ _ j^ap,upa,(Hn) ^ _ jjanupa,(&n) ^ ^ 

_ jjanvpa,(fLTi) ^ _ jja^upa, (fin) ^ _|_ £ja[i/per,(fin) ^q^m] _ £) a<T [P(())piPi' 

- l> q ^(o)p^T + # Q[CT (oKKK ] ] (9) 

i,fc,i=i 3=1 

- E s^fe - pi)"(p* - Pi) p fe - Pi) ffI Axu-i,fc-i,i-i(0) 

-E^V'CPi-Pir^SSj-iCo) 

J'=2 
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4 



+ E PiiPj -Pi) [u (Pk -Pi)V 1]a A)oj-i, fe -i(0) + pW^ 1]a ^K(0) 

j,k=2 

- E [pip^Pi - P^ p] 9 aa + PiW a fe - piT 

+ PiP l f(p 3 -Pi)V ]lQ ]^ooj-i(0) +^KK^ ]a ^oo(0). (3.17) 

As above, the coefficients -Eoojfc an d -^oooo are redundant and introduced for conve- 
nience. The coefficients of the 4th rank tensor 5-point function are then given by 

w - V^nn^M « nHfn^ 1 det(ZW) 
Arooo - - Z^^ooooW _ ^o^oooo l u ) - 48 det(F) ' 

4 

-E'oojfc = — ^ViD oo jikii'i) $ ji$ki — i] D OQ j k (0) 
i=i 

+ {sEoa^Sokw^ - ^i(o)] + (j - , 

4 
i=l 

f 4 

+ j 2 ^ Cji[-Po0fc^jmi(0^fci^^mi — A)OfcZm(O)] 
^ i=l 

+ 0'<->*0 + 0'<"> + 0' j,k,l,m= 1,...,4, (3.18) 

with 

3 3 3 3 

£>im(o)= E ^m(o) + 4 e A-w(o) + e E %(o) + 4 EA(f) + fl o(o), 

j,k,l,m=l j,k,l=l 3,k=l j=l 

3 3 3 

Am(o) = - E Dj_ 1Mm (o) - 3 E A--i,«(o) - 3 E A--i,*(o) - A-i(o), 

k,l,m=l k,l=l k=l 

3 3 

£W°) = E A--i,fc-iM( ) + 2 E4--i.*-i.'(o) + A--i,fc-i(°)> 

l,m=l 1=1 
3 

^jWl(O) = - E ^-l,fe-l,«-l,m(0) - £>j-l,fe-lj-l(0), 
m=l 

Djklm{ty = -Dj-l,ifc-l,«-l,m-l(0), 

^ffii(o) = -E^2&(o)-^ffl(o). 

-^ooooi(O) = ^ooooj-i(O)) 

3 3 3 

A)om(0) = - E ^ooifci(O) " 3 E ^ooifc(O) - 3E A»;(0) - A)o(0), 

j,k,l=l j,k=l j=l 
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3 3 

Axyn(O) = E Ax>j-i,m(0) + 2 ]T ^oo,i-i,fe(0) + Ax>j-i(0), 
fc,«=i fc=i 

3 

A)Ojfcl(O) = — ^ A)O,7-l,fc-:u(0) — -D OJ-l,fe-l(0), 

1=1 

Doo jk i(0) = Doo :j -i,k-i,i-i(0), j,k,l,m = 2,3,4. (3.19) 
4 Consistency checks, applications, and generalizations 

In order to check the explicit results for the 5-point tensor coefficients we have addi- 
tionally calculated them by applying the usual Passarino-Veltman algorithm [1], which is 
conceptually completely different from the method of this paper. In the actual compari- 
son we expressed the redundant terms Eqqq^, etc., in the Lorentz decomposition of E^, 
etc., in terms of the coefficients Ejk, etc., by exploiting the relation (3.11) for the metric 
tensor. The numerical comparison of the coefficients showed agreement between the two 
methods for non-exceptional phase-space points. 

Moreover, we have investigated the performance of the presented method in practice. 
To this end, we have implemented the new method into the calculation of the one-loop 
QCD corrections to the process gg/qq — > ttH [12] and of the one-loop electroweak cor- 
rections to e + e~ — ■> uiUiR [13], which both were originally evaluated using Passarino- 
Veltman reduction. We have observed a drastic improvement in the numerical stability 
of the 5-point function. While near the phase-space boundary the Passarino-Veltman 
approach could only be rescued by an extrapolation from the inner phase-space region 
(see Ref. [12] for details), this CPU-time-consuming procedure is practically unnecessary 
for the approach described above. The results obtained with the two methods are in 
mutual agreement within the integration errors, but owing to the more extensive use of 
the extrapolation the calculation employing the Passarino-Veltman methods takes much 
more time. 

Finally, we consider the possibility of generalizing the method of this paper to dimen- 
sional regularization. The generalization is trivial in all cases where only abelian soft 
singularities are involved. Typical examples are pure QED or electroweak processes with 
an exactly massless photon and no massless charged particles, or QCD processes that do 
not involve external gluons or massless quarks. In such cases the soft singularity arises 
from diagrams with photon or gluon exchange between two external (on-shell) massive 
lines and shows up as a single pole in e = (4 — D)/2. Alternatively, if for D = 4 an 
infinitesimal photon or gluon mass A is chosen as IR regulator, the singularity leads to 
In A terms. The correspondence between these regularizations is well known (see e.g. 
Ref. [11]): 

mi * — * (wa+e) +0(e) (41) 

If external massless charged particles are involved but no external massless gluons, 
i.e., if the IR singularities are due to photons or gluons coupled to massless fermions 
or sfermions the following approach can be used. The whole calculation can be carried 
out with mass regulators obeying the hierarchy A <C m <C \Q\, where A is again an 
infinitesimal photon or gluon mass, m is a small fermion or sfermion mass, and Q denotes 
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a typical mass scale of the process. Then these results can be translated into dimensional 
regularization as described in Ref. [11] for the complete QCD and SUSY-QCD amplitudes. 

If non-abelian soft singularities or overlapping soft / collinear singularities are involved 
the situation is more complicated. A convenient possibility to make use of the 4-dimen- 
sional approach of this paper is, for instance, described in Ref. [12]. There, a method 
is presented for translating D-dimensional into 4-dimensional integrals by constructing 
regularization-scheme-independent finite integrals upon subtracting well-defined, simple 
auxiliary integrals with the same singularity structure. These auxiliary integrals are 
entirely built from 3-point functions. In summary, this means that D-dimensional 5-point 
integrals are first converted into 4-dimensional 5-point integrals (regularized with masses 
or off-shell momenta) and 3-point integrals. The 4-dimensional 5-point integrals are then 
decomposed into 4-point integrals with the method described in this paper. 



5 Conclusion 



A method for reducing one-loop tensor 5-point integrals to related standard 4-point 
integrals is proposed that entirely avoids inverse Gram matrices, which are potential 
sources of numerical instabilities in practice. The presented explicit results for tensor 
coefficients of 1-, 2-, 3-, 4-, and 5-point functions comprise all expressions needed to deal 
with 1 — * 2, 2 — > 2, and 2^3 particle reactions with up to four external gauge bosons of 
renormalizable gauge theories and all expressions, apart from 6-point functions, for 2^4 
particle reactions with up to two external bosons. The relevant 5-point functions are UV 
finite, but may contain IR (soft or collinear) singularities. In the explicitly given results a 
four-dimensional regularization scheme is assumed, but possible ways of translating them 
into dimensional regularization are described. 



Appendix 

A Calculation of extra terms resulting from UV-divergent 4-point functions 



For the reduction of E pvp and E pvpa we need to evaluate the UV-divergent 4-point 
functions D^ up and D^ upiT introduced in Section 2. Recall that these integrals are defined 
in four space-time dimensions. Going to D dimensions and using the large-mass expansion 
[15] for the regularization parameter A, we obtain for these integrals 



apvp 



(2tt/x) 



4-D 



17T 



d D q 



-A 2 



AyVxAWa q 2 - A 2 



D=4 



(27T/X) 



T S~ I dDq q tr P ? - A 2 + 0(A ' 2 ln A) + ° {D ~ 4) (A - 1} 



and 



D (A) 
apvpa 



(2tt/x) 



4-D 



d D Q a qpquq P qa 



-A 2 



A^oo Dap,vpa + 



NqNiN 2 N 3 q 2 - A 2 
(27r/i) 4 - D f^ D qgq^qpq^ 



17T Z 



D=4 

-A 2 



(g 2 ) 4 q 2 - A 2 



1 iqpi 

2^ 2 

z=o q 



+ 0(A- 2 lnA) +0(D- 4), 



(A.2) 
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where D afiup and D afiupa are the usual 4-point functions in dimensional regularization as 
defined in (2.1), and we included a non- vanishing momentum p in N Q for later conve- 
nience. The extra vacuum integrals can be evaluated by standard methods; the results 
are 



(27T/i) 



4-D 



(2tt//) 



4-D 



17T 



/ 

/ dD < 



-A 2 



(g 2 ) 4 g 2 -A 2 
-A 2 



(q 



2\4 



q 



A 2 



1 3 

If 



96 ^ 



«=0 



2^9a[p9vp] 

1 y 2 ^ 

z=o y 

2 



A + ¥ + ln A^ 



A+ 12 +ln A^ 



+ 0(D-4), 



(A.3) 



where 



7e + In Air 



(A.4) 



4-D 

represents the UV singularities in dimensional regularization, 7 E is Euler's constant, and 
\x is the dimensionful parameter of dimensional regularization. 

Since the UV singularities of D^ Mp are regularized by A, the poles in D — 4 in the 
vacuum integrals cancel the poles in D apup and D apiUpa , and we can write 



L> (A) 



£)(fin) + A 

fix. .M P A^oo m-vp 



ii\...v P 



(A.5) 



with 



A 



pupa 



A 



apvpa 



1 

24 



9n[v9pcr] 



11 , M 2 



+ 0(D-4), 



1 C 

96 v 



X^Pi)^^^^]] 



z=o 



25 , /i 2 



+ £>(£- 4) 



(A.6) 



and 



n(fin) 

pupa 



n(fin) 

apupa 



D pv( „ - —g^g^A + 0(D-4), 
1 3 

D apupa + — — ^) [a9p[u 

g pa]] A + 0(D-A). 



(A.7) 



The quantities Up , D^ Mp , and A Ml ...^ p are finite in £) dimensions, so that we can 
perform the limit -D — > 4 there. 

In the reduction of -E^p, UV-divergent contributions appear only in the difference 
[A apup {i) — A apup (0)]. Since A apiVp is momentum independent, the difference, and thus 



Ufj, up , vanishes. 



For E^pty, we have contributions from A MI/p<T (0) and A ailV( „{ 
terms are given by 



'■apupa 



(0). These 



11 , /i 2 
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1 4 

^anupc (i) - ^apupa (0) = — (J^ Pi) [a9p[u 



9pa]} 



1=1 



25 /i 2 
12 +ln A^ 



1 4 



25 , /i 2 



ggP«,[o;fl l //[^5'pCT]] 

Inserting (A. 8) into (2.14) and using 



25 /x 2 ' 
12 +ln A^ 



+ £>(£> -4) 
^O(D-A). (A.8) 



E det(z£>) 2p? ft>a = E det(Z( 4 )) (zW)r.i 2p? Pi , a = 4det(Z( 4 )), (A.9) 

which employs the relation between det(zff) and the inverse Gram matrix (Z^) -1 , we 
find 



£W = -^^]<rdet(ZW) + 0(jD _ 4 ). 

B Reduction of tensor 1-, 2-, 3-, and 4-point functions 

One-loop tensor iV-point integrals have the general form 



(A.10) 



T^.-.ppiPu -,Pn-i, m , m^v-i) 



(2tt/x) 



4-D 



m 2 



J q N N 1 



N N .! 



(b.i; 



with iVj, i = 0, ...,N — 1, defined in (2.2). Following the notation of Ref. [9], i.e., 
T 1 — > A, T 2 — > _B, T 3 — > C, T 4 — > and using the conventions of Ref. [5], we decompose 
the genuine tensor integrals into Lorentz-covariant structures: 

A» = 0, = g^A 00 , 

B"=^B U = g^B 00 +p^ l B ll , 

B^ up = 9^ u PiB 001 + PiPiPiB ul , 



c» = Y.PjC 3 , 

3=1 
2 



j,k=i 



c^ = Y.9^pfc mj + E tfrWCw, 



3=1 



j,k,l=l 
2 



cm = g^gd'Cum + E (9^p?pI + g^p^Coojk + E tirtrftinC jiam , 

j,k=l j,k,l,m=l 
3 3 

D» = E Ppj, D*" = gTDw + E PjP U k D 3k, 

3=1 j,k=l 

= j^^pfDooj + E P^Djki, 

3=1 3,k,l=l 
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jywp* = g^gP]°D 0000 + £ (^Pfpl + g^p-pi^k + E P»fp a m D jk lm, 

j,k=l j,/c,J,m=l 
3 3 

= Y,g [[llp g p]a p T ^D om , 3 + E (^p&>[ ] +^^r I )^«wH 

J=l j,k,l,=l 

+ E P^p a mP T nD jklmn . (B.2) 
j',fc,i,m,n=l 

Because of the symmetry of the tensor all coefficients Bi,...,DjM mn are symmetric 

under permutation of all indices. They are reduced to scalar integrals recursively by 
inversion of systems of linear equations. The inhomogeneity of these equations consists 
of coefficients of lower rank [1]. The equations of this system are obtained by contracting 
one integration momentum with the (N — 1) external momenta pf 1 and for N > 2 
also by contraction with the metric g^ lfJ/2 . Using 

2 Pi q = Nt - N - fi with f i =p 2 i -m 2 i +m%, (B.3) 

in (B.l), the first two terms on the right-hand side of (B.3) each cancel exactly one 
propagator denominator of q^T^ p , the third term is proportional to plp . Likewise 
the contraction with g^ 1 ^ 2 yields a factor q 2 in the numerator of g^^T^ , which can 
be written as q 2 = N Q + m 2 ,. The N term cancels the first propagator, the second term 
leads to the tensor T^ 3 p . With the abbreviation T^ 1 (i), denoting p with the 
ith denominator omitted, this yields 



UlrpN =T N-1 (j\_rpN-l (n\ _ f.rfN 



g^K^ p = T^ p (0) + m 2 T^ p . (B.4) 

Note that for T^' 1 (0) a shift of the integration momentum q^ — * q^ —p^ has to be done 
in order to achieve the standard form (B.l). The tensor integrals with shifted momenta 
T / ^~J tp (0) are defined as in (3.4). Expressing the tensors in (B.4) by (B.2), the desired 
recurrence relations can be read off by comparing coefficients. 

In the following we summarize the results for the reduction of all tensor coefficients 
defined in (B.2). These results can also be read off from the generic results in Section 4.2 
of Ref. [5]. Note that several tensor coefficients can be obtained from different reduction 
formulas. This allows to check the correctness and the numerical stability of the results. 

The tensor coefficient of the relevant 1-point function reads 

A 00 = ^m 2 A + Knl. (B.5) 



For the tensor coefficients of 2-point functions we find 

1 

2p{ 
1 



B l = —[A (l)-A (0)-f 1 B ], (B.6) 



-Boo 



6 

1 

3p? 



A)(0) + fxBx + 2m 2 B + m 2 + m{- \p\\ , 
A)(0) - 2j x B x - m 2 B - \ (m 2 + mj - , (B.7) 



17 



Am = 7TT [^oo(i) - A)O(O) - hB t 



2p 



ooj 



B 



2mgBi - A o (0) + /i5 n - ±(2mg + 4m? - p\) 
1 



111 



4pf 



A (0) + 3/iSn + 2mgS! - | (2m^ + 4m? - pfj 



For the 3-point functions we obtain for the vector case 



c = E(z (2) ),"»«i. -1.2. z «=(^fM 



n=l 



with 



< = S (n) - B (0) - f n C , n = 1, 2, 
for the second-rank tensor case 



Coo = |^oCo + i-Bo(O) + | X! /m^™ + I' 

m=l 

C7 y = E(Z< 2 >)ri[i$. - 2C oM, J = 1, 2, 

n=l 

with [i„ is defined in (3.5)] 

R 2 ni = B 1 (n)5 nt - B t (0) - f n C h n,i = l, 2, 

and 

Bx(0) = -Sx(0) - S o (0), fl 2 (0) = 5x(0), 
for the third-rank tensor case 

2 

Cooi = ImlCi + |-Bj(0) + | X! fj^ij ~ T8> 

3=1 

2 

C 0i = ^2( Z ^)in R n00i 
n=l 

C ijk = J2( zi2) )Tn[Rljk - ZCooAk - 2CW<y , k — 1,2, 



n=l 



with 



-^noo = B 00 (n) - B 00 (0) - f n C 00 , 

R lii = B n (n)5 ni 5 nj - Bij(0) - f n C ij: n, i,j = 1, 2, 
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and 



B u (0) = 5ii(0) + 2B 1 (0) + B o (0), 
B 12 (0) = -BuW-MO), 
B 22 (0) = B n (0), 
and for the fourth-rank tensor case 

2 

Coooo = jm 2 C o + |-Boo(0) + § ^ fiC 00 i + js( m o + m l + m l) ~ m [Pi + #2 + (Pi 



i=l 
2 



fe=i 



n=l 
2 



CijfcZ — X]('^^)m 1 [-^njfei ~~ 2C 0jk8nl ~ 2C 0fc^nj ~ 2C o/j5nfc], i,j, k, I — 1, 2, 

n=l 

with 

-^nOOi = -SoOl(^)^ni _ -BoOi(O) — / n CoOi, 

-^nijfc = ^lll(^)4i4j^nfe --Bijfc(O) ~ fnQjk, n,i,j,k = 1,2, 



and 



Sooi(O) = -Sooi(O) - 500(0), 5oo2(0) = S OO i(0), 

£ m (0) = -S m (0) - 3S U (0) - 3^(0) - B (0), 
B 112 (0) = B m (0) + 25 n (0) + Si(0), 
#122(0) = -B lu (0)-B n (0), 
B 222 (0) = S m (0). 

The tensor coefficients of 4-point functions are given by 

3 / 2pipi 2pip 2 2pip 3 N 

A = E( Z(S) )M, i = 1,2,3, ^ 3 ) = 2p 2Pl 2p 2 p 2 2p 2 p 3 

\2p3Pi 2p 3 p 2 2p 3 p 3 / 



n=l 



with 

Sn = C (n) - C o (0) - /„£>„, n = 1, 2, 3, 

for the vector case, by 

3 

Dqo = m 2 D + iC (0) + ± X / m An, 

m=l 



Aj = E(^ (3) )m[^ - 2£>oo^], i, j = 1, 2, 3, 
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n=l 



with 

S 2 m = C in (n)$m ~ C^Q) - f n Di, n,i = l, 2, 3, (B.23) 

and 



C 1 (0) = -X;C i (0)-Co(0), a(0) = C^i(0), * = 2,3, (B.24) 

i=l 

for the second-rank tensor case, by 

3 

D 00l = \rn\D, + \Cm + \ fjDij, 



j'=i 



a™ = E(^ (3) )r^oo, 

n=l 



A, fc = E(^ (3) )r„ 1 [^-2 J Doo^n fc -2 J D 00fc M, i,j,fc = 1,2,3, (B.25) 



n=l 



with 



^00 = CooW- Coo (0) "/„£>, 



00) 



= C inj „(n)(J ni (J nj - C^(0) - / n Ai, n, i,j = 1, 2, 3, (B.26) 



and 



2 2 

C n (0) = J2 c ij (o) + 2jrc i (o) + c (p), 

i,j=l i=l 

Cii(o) = -E^(o)-^i(o), 

J'=l 

C,,(0) = A-i,,_i(0), i,j =2,3, (B.27) 
for the third-rank tensor case, by 

3 

-DoOOO — 3^0 A)0 + 6^0o(0) + gE /iA)0i + 3g, 

i=i 

3 

-Dooij = \ m oDij + |Cy(0) + \ fkDijk, 

k=l 

3 

A)0ij = XK^^i^f^nOOj _ 2D 000^nj], 



n=l 



Ajm — y^(-^^)m 1 ['5'n-7-fc/ ~~ 2D 00 j k 5 n i — 2D om 5 n j — 2D oij5 n fe]) J, k,l — 1, 2, 3, 

(B.28) 



n=l 
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with 

SnOOi = Cooi n (n)5 n i — Cooi(O) — f n Dooi, 

s nijk = C injnkn (n)5 ni 5 nj 5 nk - (7^(0) - f n D ijk , n, k = 1,2, 3, (B.29) 

and 

Cooi(O) = ^00,(0) - Coo(0), C 00t (0) = C O o,*-i(0), 
i=i 

Cm(0) = - E Ciifc(O) - 3 E Cii(O) -3^,(0) - C o (0), 

i,j,k=l i,j=l 1=1 

2 2 

c ni (o) = E A-i Jfc (o) + 2Ea-i J (o) + A-i(o), 

j,fc=i i=i 

2 

Ciij(o) = - E Ci-i,j-i,fc(o) - Ci_ij_i(o), 
fe=i 

Cy fc (0) = Ci-u-Lk-^O), i, = 2,3, (B.30) 
for the fourth-rank tensor case, and by 

3 

-DoOOOi — 4 m 0-^00i + gCoOi(O) + g E fj^OOij ~ 192) 



Dooijk — \w?§Dij k + |Cjjfc(0) + § E fmDijkm, 

m=l 

3 

-DoOOOi = E(^ (3) )^ 1,S 'n0000> 
n=l 

3 

-DoOijfc = E(^ (3) )m 1 [^n00jA: ~ 2D 000j^nfc ~ 2D 0000k 5 nj ] , 



n=l 
3 



Dijklm — ^2(Z^)in[Snjklm ~~ ^^OOjkl^nm ~ 2D 00klm 5 n j — 2D 00 i m j5 nk — 2D 00m j k 5 n i], 
n=l 

k, I, m — 1, 2, 3, (B.31) 

with 

^noooo = Coooo(^) — Coooo(O) — f n D 0000 , 
Snooij = Cooi n j n (n)5 ni 5 n j — C O oy(0) — /n-Dooij, 

^nijfcz = C in j nkn i n (n)5 ni 5 n j5 nk 5 n i — (7^(0) — f n D ijk i, n, i,j, k, I — 1, 2, 3, (B.32) 

and 

2 2 

Coon(O) = E ^(0) + 2^^00,(0) + C 00 (0), 
«j'=i i=i 
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CoOil(O) — — ^Coo,i-lj(0) — CoO,i-l(0), 
3=1 

Cooij(O) = Coo,j-ij-i(0), 

2 2 2 2 

Cmi(0)= E C m (0)+A J2 ^(0)+6^4(0)+4^Q(0)+C o (0), 

i,j,k,l=l i,j,k=l *,i=l i=l 

C mi (0) = - f; Q-ijfcKO) - 3 2 Ci-ijfc(O) -3^(7^(0) - Ci-i(O), 
j,fe,/=i j,fc=i j=i 

2 2 

di«(o) = E a-ij-i,«(o) + 2x;a-ij-i,fc(o) + c' i _i J -_i(o), 



fc,«=i 

2 



k=l 



Cujk(0) — — J~] Cj-ij-i,k-u(0) — Cj-i,j-i,fc-i(0), 

Cijfci(O) = Ci-i,j-i,ifc-i,i-i(0), i,j,k,l = 2,3, 
for the fifth-rank tensor case. 



(B.33) 



C UV-divergent parts of tensor integrals 

For practical calculations it is useful to know the UV-divergent parts of the tensor 
integrals explicitly. We directly give the products of D — 4 with all divergent one-loop 
tensor coefficient integrals appearing in renormalizable theories up to terms of the order 
O(D-A) 



(D-4) A (m 
{D - 4) A 00 (m 
(D - 4) 5 (pi,m ,mi 
(D - 4) B 1 (p 1 ,m ,m 1 
(D - 4) B 00 (p 1 ,m ,m 1 
(D - 4) Bn(pi,ffl ,mi 
(D - 4) B 00 i(pi,m ,m 1 
(D - 4) £m(pi,m ,TOi 
(D - 4) C 00 (pi, P2, mo, mi, m 2 
(D - 4) C oi(pi, P2, m Q , mi, m 2 
(D - 4) Coooo(pi, P2, m , m 1 , m 2 

(D - 4) Cooij(pi,P2,mo,mi,m 2 
(D - 4) Doooo(pi,P2,P3,mo,m 1 ,m 2 ,m 3 



-2ml, 



-\rriQ, 
-2, 
1, 

\(pl - 3ml - 3ml), 

2 
"3' 

-i(p?-2r^-4m?), 



i 

2' 
1 

'2' 
1 

6 ■ 



- ^{ml + ml + ml), 



j_ 

12' 
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(D — 4) D Q0Q0i (p 1 ,p 2 ,P3, m ,mi,m 2 , m 3 ) = ^. 
All other scalar coefficients denned in (B.2) are UV finite. 



(C.l) 
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